NEW RESULTS ON TOPOLOGICAL SUSCEPTIBILITY IN SU(3) GAUGE THEORY 
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We survey recent lattice results on QCD topological properties. The behaviour of the topological susceptibility 
at the deconfining phase transition has been determined. This advance has been made possible by an i) an 
improvement of the topological charge operator and ii) a non-perturbative determination of renormalizations. 
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1 Introduction 

Relevant progress in the study of topological prop- 
erties of QCD has been made possible by two re- 
cent developments in lattice gauge theories: 

1) Improvement of the operatorscl. 

2) Non pertnrbative determination of the renor- 
malizations □. 

To clarify the meaning of 1) we recall that the 
building block of lattice gauge theories is the link, 
Uu{n) — exp(zaA M (n)), or parallel transport from 
the site n of discretized spacetime to the neigh- 
bouring site in the direction (X, ft + aft (a is the 
lattice spacing). All the operators are constructed 
in terms of links, and therefore contain arbitrar- 
ily high powers of aA^. The identification with 
continuum is usually done by requiring that the 
leading order in the power expansion of the lattice 
operator Ol concides with its continuum counter- 
part O 

L = Oa d + 0(a d+1 ) (1) 

d is the dimension of the operator. Higher order 
terms in eq. ([!]) can be changed with large arbi- 
trariness. 

For the action, for example, eq. ([[]) reads 

S L = -^G (lu G^a i + 0(a 6 ) (2) 

Wilson's action = £(1 - 1/3 Re TrD) obeys 
eq. (§). 

The idea of improvement consists in exploiting 
the arbitrariness in higher order terms to reduce 
lattice artifacts. Improving the action can make 
the lattice larger in physical unitsQ. Our approach 
will be to keep the usual Wilson action and to i 
prove the operator Ql for the topological chargel 
to reduce renormalizations from lattice to contin- 
uum. Renornijh^tions will be determined non 
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2 Topology in QCD. 

The key relation is the anomaly of the singlet axial 
current 

= 2N f Q(x) (3) 

where 

fi = J2^fw 5 Tpf ( 4 ) 

/ 

/ means flavour, and Nf number of flavours. Q(x) 
is the topological charge density 



Q{x) = 



JL_e^G° v (x)G a pa (x). (5) 



The anomaly can explain the magnitude of 
the rj' mass (the so called U{1) problem^) if the 
topological susceptibility of the vacuum in the 
quenched approximation (leading order in ■A- ex- 
pansion) 



\ = / d i x{0\T(Q(x)Q(0))\0). 

is large enough. Quantitatively Bi 

2#/. 
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(6) 



(7) 



Testing eq. (0) is a check at the same time of 
QCD and of the -A- expansion, and can be done 
on the lattice. 

The behaviour of x at the deconfining temper- 
ature is .also an important test of models of QCD 
vacuum 13. 

Another use of eq. ([?]) which can be made on 
the lattice is the measurement of the singlet axia 
charge of the nucleon, Gi(0), which is given bytL2 



Gi(0) ~ lira 



Nf tfs'\Q\ps) 
fXff M p u s >(p')j 5 u s (p) 



(8) 



perturbativelyt 



G\ is usually related to the quark spin content of 
the proton til. 
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3 Topology on the lattice. 

A lattice version of the topological charge density 
operator is 



Xx) 



-1 ±4 

^2 



29^2 L^i 



(9) 

In the formal limit a — > 0, in accordance with eq. 
(1 

Q L {x)^a i Q(x) + 0{a 6 ). (10) 

Any other choice for Ql will differ from (j^) by 
0(a e ). 

Ql{x) renormalizcs multiplicatively 

Q L (x) ~ ^ x Z((3)Q(x)a(f3) 4 + 0(a e ) (11) 

Similarly, for the lattice topological suscepti- 
bility, defined as 



Xl = (]TQl(x)Ql(0)), 



(12) 



we have 







XL 

M{fi) 



p^ZipfaiPfx + Mifi) (13) 
B(f3)(G 2 )+P(P){l)+(D(a e ) 



where G2 is the trace of the energy-momentum 
tensor and 1 is the identity operator. The pres- 
ence of M(j3) is due to the fact that the lattice 
regularization does not obey the prescription for 
the..&ingularity at x = which defines x of eq. 
0)DH. From eq. ([[J) 



X 



XL - M(J3) 
Z 2 a A ■ 



(14) 



Xl is determined numerically on the lattice. 
P{/3), Z{0) and B{0) depend on the choice of Ql, 
i.e. on the specification of the terms 0(a 6 ) in eq. 
(|Io|). A good choice gives a small M((3) and Z w 1, 
so that in eq. (^4|) x~Xl and lattice artifacts are 
a small part of the numerical determination. 

To determine Z, Ql = J <\ a xQl(x) = ZQ is 
measured on an ensemble of configurations with 
a definite value of Q. A one instanton configura- 
tion is heated to dress it with short range quantum 
fluctuationSj-Thc number of instantons is checked 
step by stepE3, and the configurations in which ei- 
ther the initial instanton has disappeared or new 
instantons have been created are discarded from 



the sample (Fig. 1). M((3) is determined by mea- 
suring xl on an ensemble of configurations with 
Q = (0 instantons plus quantum fluctuations). 
Here again the global topological charge is checked 
step by step and configurations in which it has 
been changed in the heating procedure are dis- 
carded (Fig. 2). We have used an improved oper- 
ator obtained from Ql(x) by successive smearings, 
and which differs from it by terms C(a 6 )El. 

Figure 3 shows the plateaux reached by heat- 
ing 1 instanton and the Z's for the 0, 1 and 2 
improved operators. 




Figure 1: Distribution of topological charge Q for a set of 
2000 configurations obtained by 15 heat-bath updatings of 
a 1-instanton configuration, fi = 5.75. 




Figure 2: Distribution of topological charge Q for a set of 
500 configurations obtained by 36 heat-bath updatings of 
the zero field configuration, fi = 5.90. 



The analogous behaviour for M{0) is shown in 
Figure 4. Improving twice the operator produces 
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a factor of 3 in Z, or a factor of 10 in Z 2 , and a 
factor of 10 reduction of M((3). The ratio between 
the physical signal and the lattice artifacts is then 
improved by 2 orders of magnitude going from the 
to the 2 smeared operator, and the subtraction 
in eq. ( |l4|) becomes of the order of 10% of \l in 
the scaling window. 
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Figure 5: Determination of x at T = 0. Diamonds, squares 
and circles correspond to 0, 1 and 2 smearings respectively. 



Figure 6 shows a new result, which has been 
made possible by the improvement: the behaviour 
of x across T c . The \ drops by at least one order 
of magnitude at T c . This result can provide an 
important check of models of QCD vacuumu. 



Figure 3: Determination of Z by heating of a 1 instanton 
configuration. Diamonds, squares and circles correspond 
to 0, 1 and 2 smearings respectively, f) = 6.36. 
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Figure 4: Determination of M(/3) by heating of the flat 
configuration. Diamonds, squares and circles correspond 
to 0, 1 and 2 smearings respectively, f) = 6.30. 



Figure 5 shows the result for x at T = deter- 
mined with the 0, 1 and 2 smeared operators: they 
agree with each other, scale properly and confirm 
previous results in favour of the Witten-Veneziano 
mechanism crli3. 




Figure 6: Determination of \ at T ^ 0. Circles and squares 
correspond to 1 and 2 smearings respectively. The horizon- 
tal line is the T = result. 



Difficulties are encountered in thermalizing the 
topological charge with the hybrid Monte Carlo al- 
gorithm in full QCDEJ (Figs. 7,8,9): we are, how- 
ever, implementing a multicanonical algorithm in 
which the quark mass is used as a new variable 
in the Monte Carlo, to exploit the fact that for 
high quark masses this difficulty is less severe. The 
procedure works, so that the same technique de- 
scribed here will soon allow the determination of 
X, x' in full QCD and the measurement of the spin 
content of the proton. 
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Figure 7: Time history, in units of molecular dynamics 
time, of the topological charge Q for a hybrid Monte Carlo 
simulation at /3 = 5.35 and am = 0.01. 
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Figure 8: Distribution of values of Q in a sample of 1000 
configurations in full QCD obtained with the hybrid Monte 
Carlo algorithm. (Q) ^ and the distribution is not sym- 
metric under Q — > — Q in contrast with what happens with 
usual heat bath algorithms for quenched QCD, see Fig. 9. 



Figure 9: Distribution of values of Q in quenched QCD 
obtained with the heat bath algorithm. 
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